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Abstract. We use recent population synthesis results to investigate the distribution of pulsars in the frequency 
space, having a gravitational strain high enough to be detected by the future generations of laser beam interfer- 
ometers. We find that until detectors become able to recover the entire population, the frequency distribution 
of the 'detectable' population will be very dependent on the detector noise curve. Assuming a mean equatorial 
deformation e = 10~^, the optimal frequency is around 450 Hz for interferometers of the first generation (LIGO or 
VIRGO) and shifts toward 85 Hz for advanced detectors. An interesting result for future detection stategies is the 
significant narrowing of the distribution when improving the sensitivity: with an advanced detector, it is possible 
to have 90% of detection probability while exploring less than 20% of the parameter space (7.5% in the case of 
e = 10~^). In addition, we show that in most cases the spindown of 'detectable' pulsars represents a period shift 
of less than a tens of nanoseconds after one year of observation, making them easier to follow in the frequency 
space. 
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1. Introduction 

The coming years will be marked by the beginning of op- 
erations for the major gravitational wave interferometric 
antennas LIGO, VIRGO, GEO and TAMA. Starting with 
a sensitivity of the order of h w 10^^^, they are expected 
to evolve in a few years into second generation experi- 
ments with a sensitivity improved by one or two orders of 
magnitude depending on the frequency. Ground based de- 
tectors will scan the sky searching for gravitational waves 
with frequencies between tens of hertz up to tens of kilo- 
hertz. Potential sources fall roughly into three classes: 
bursts, stochastic background, continuous, involving dif- 
ferent search techniques as match filtering for coalescing 
binary systems and cross-correlation between detectors for 
the stochastic background. The detection of pulsars (and 
other nearly monochromatic continuous sources) can be 
achieved by integrating the signal during times of about 
10^ s and searching for statistically significant peaks at 
fixed frequencies in the power spectrum. This method be- 
comes quite complicated as soon as one consider that: 
1) pulsars are not really monochromatic emitters. The 
Doppler effect due to the Earth motion and the intrinsic 
spindown of the star become significant over long period 
of times spreading power across many frequency bins and 
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must be corrected before performing the fourier transform. 
2) most of pulsars are not observed in the electromagnetic 
domain due to strong selection effects, making necessary 
an all-sky and all-frequency search. The computational 
cost is therefore the main concern of searching for contin- 
uous sources: if the whole parameter space is scanned, the 
computing power required to process the data will be out 
of any reasonable value: of the order of 3 x 10^" Tera-flops 
for integration times of one year, frequency bandwiths of 
2500Hz and decay times of 1000 year (Frasca et al. 2000). 
Alternative search algorithms involving hierarchical meth- 
ods which follow up candidate detections from a first pass 
search have been investigated during the past few years 
(Brady et al. 2000, Schutz & Papa 1999, among others). 
These techniques reducing very significantly the comput- 
ing time, by more than ten orders of magnitude for a small 
loss of sensitivity (by less than an order of magnitude) 
(Frasca et al. 2000) represent the best hope of detection 
for actual computational ressources and detector sensitiv- 
ities. An attempt to optimize these algorithms could be 
the determination of the most likely range on the parame- 
ter space. Without restricting the search area definitively, 
the probability to have a first detection early could be op- 
timized. Brady et al. (1998) have considered a survey in a 
limited sky area, searching for pulsars in the galactic core. 
In a previous work (Regimbau & de Freitas Pacheco 2000, 
hereafter paper I), population synthesis of galactic normal 
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pulsars was performed in order to recover the statistical 
properties of the real population. This approach permits 
the modeling of the galactic pulsar contribution to the 
gravitational emission and to deduce its detectability by 
laser beam interferometers. The gravitational strain being 
proportional to the star equatorial deformation, the num- 
ber of detection depends on the mean value adopted for 
the equatorial ellipticity. This parameter is quite uncer- 
tain since it depends on the cooling history and on differ- 
ent strains underwent by the crust. Past studies (Thorne 
1980) have shown that equatorial distortions cannot be 
larger than 10~^, expected to be an upper limit related 
to the maximum stress that the solid crust could support, 
but it is worth mentioning that it could be several or- 
ders of magnitude smaller. Our simulations have showed 
that 10~^ is the critical value to have at least one de- 
tection with interferometers of the first generation (10~^ 
for an advanced interferometer) In this paper, by using 
our previous population synthetis we investigate the de- 
tection probability of a pulsar in the frequency space. In 
order to illustrate our findings, we report here results for 
mean ellipticities equal to 10~^ and 10~^, as well as how 
the results are affected by the detector sensitivity. 



2. The Population Synthesis 

2.1. The Model 

Population synthesis aimed to estimate the gravitational 
contribution of galactic normal pulsars are described in 
details in paper I. The initial parameter distributions be- 
ing fully taken into account, our statistical approach gives 
more realistic estimates of the gravitational strain and the 
number of 'potentially detectable' pulsars with frequen- 
cies falling into the detector sensitivity band, than other 
analytical methods. Among the main results of our simu- 
lations those rclevants for the present purpose are briefly 
reminded here: a) the average initial period of pulsars is 
290 ms with a dispersion of 100 ms. Pulsars are not born in 
majority as fast rotators in agreement with recent studies 
(Bhattacharya ct al. 1992, Hartman ct al. 1997, Narayan 
1987) and with observations, since the mean period of pul- 
sars younger than 20, 50 and 100 kyr is respectively 123, 
162 and 173 ms (Lyne & Granham-Smith 1998). However, 
the present simulations with a non-fixed initial period per- 
mit to account for the fastest objects like the Crab pulsar 
(33.4 ms) and for the diversity of young pulsar periods: 
PSR B1509-58 (as young as the Crab pulsar) has a pe- 
riod of 150 ms, PSR B1610-50 (six times older) a period 
of 231 ms, while the VELA pulsar and PSR B195H-32, 
respectively 10 and 100 times older than the Crab pul- 
sar rotate with periods of 89 ms and 59 ms. b) around 
150000 pulsars are expected to be present in the Galaxy, 
c) among them, pulsars in the detector frequency range of 
sensitivity (called later "potentially detectables" ) would 
be around 1000 if the lower bound is 10 Hz or around 
5200 if it is 5 Hz. The latter is expected to be attained 
in a second phase of the VIRGO experiment, thanks to 



a sophisticated pendulum system called superatenuator. 
Notice that the upper bound of the detector (around 5-10 
kHz) doesn't interpose on the detection of pulsars, since 
the Keplerian stability limit of neutron stars corresponds 
to gravitational wave frequencies of 2-4 kHz, depending 
on the equation of state. 

2.2. The Detectable Population 

When any precessional effect is neglected the two polar- 
ization components of gravitational waves emitted by a ro- 
tating neutron star are (Bonazzola & Gourgoulhon 1996) 



h+{t) = 2A{1 + cos^ i) cos(20t) 
and 

h^{t) = 4A cos i sm{2flt) 



(1) 



(2) 



where i is the angle between the spin axis and the wave 

propagation vector, assumed to coincide with the line of 
sight, and the wave amplitude is defined by 

G 



rc-* 



(3) 



In this expression, G is the gravitation constant, c is the 

velocity of light, r is the distance to the source, fl is the 
angular rotation velocity of the pulsar and the ellipticity 
£ is defined as 



e = 



1-xx lyy 



(4) 



The gravitational strain amplitude induced in the detector 
depends also on the direction and polarization of the wave 
with relative to the detector arms and it is given by: 



h{t) = h+{t)F+ {e, cj>,i,) + h^ {t)F^ {e, cj>, ij) 



(5) 



where f+.x are the beam factors of the interferometer 
(Jaranowski ct al. 1998), which are functions of the zenith 
distance 0, the azimuth (j) as well as of the wave polariza- 
tion plane orientation ip. 

Among the "potentially detectable" population men- 
tioned above (pulsars with periods less than 0.4 s for 
VIRGO and less than 0.2 s for LIGO), 'detectable' objects 
are pulsars with gravitational strain amplitudes above the 
antenna sensitivity curve. The detectability criterium can 
be expressed as follows: 



hnjfg. 

Vt 



(6) 



where h is the gravitational strain, T the integration time, 
hn the detector sensitivity (in Hz~^/^). In our analysis, 
the planned LIGO II was used as an exemple of second 
generation detectors. The factor two corresponds to the 
adopted signal-to-noise ratio. As one should expect, the 
less is e, the harder is the detectability condition. The 
number of detections as a function of the ellipticity and 
detector sensitivity is reported in Table 1. An average el- 
lipticity of 10~^ represents the detectability threshold for 
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Table 1. Number of detections as a function of the el- 
lipticity for interferometers of the first generation (LIGO 
and VIRGO) and for an advanced interferometer 



e 


10"^ 


10"® 


10-^ 


VIRGO 


15 


3 





LIGO I 


12 


2 





Advanced 


90 


12 


2 



detectors of the first generation as VIRGO or LIGO I: 
around 12—15 detections are expected if£ = 10~^, 2 — 3 
detections if e = 10~^ and no detection if e < 10~^. For 
an advanced interferometer this threshold corresponds to 
a mean cllipticity of 10^^. This result has quite interesting 
consequences: even in the case of no detection, an upper 
limit can be set to the cllipticity and improvements in the 
sensibility will gradually push down such a limit. 

As shown in Table 1, detectable pulsars are expected 
to be few at the beginning, since only the nearest and the 
fastest objects have strain amplitudes high enough to be 
above the noise level. The detector sensitivity acts as a fil- 
ter affecting a priori both spatial and frequency statistics. 
Actually our simulations have shown that the distance 
distribution is quite similar for the detectable population 
and the total population, with a peak around the galac- 
tic center. As a result the sky direction distributions do 
not depend significantly on the antenna sensitivity or on 
the ellipticity. On the opposite, the frequency distribution 
is very dependent on both of them and requires a more 
detailed analysis. 

2.3. The Gravitationnal Frequency 

The gravitational frequency distribution of detectable pul- 
sars as a function of the antenna sensitivity, for two dif- 
ferent values of the ellipticity, the theoretical upper value 
e = 10~^ and the detectability limit for interferometers of 
the first generation, e = 10~^ are shown respectively in 
Fig. 1 and Fig. 2. 

The sensitivity improvement enables the detection of 
slower objects, so the maximum of the detectable fre- 
quency distribution shifts toward lower values. For e = 
10~^, the maximum is respectively around 180 Hz and 
135 Hz for VIRGO and LIGO I and decreases to 85 Hz 
for an advanced interferometer. For e = 10~^ it is 480 
Hz for VIRGO, 435 Hz for LIGO I and 85 Hz for for an 
advanced interferometer. The difference between the opti- 
mal frequencies of the two detectors of the first generation 
arises from the fact that LIGO I is more sensitive than 
VIRGO between 50 Hz and 330 Hz and less sensitive else- 
where. What is interesting to notice for future detection 
stategies, is the significant reduction of the distribution 
width between the first and second generation (see Table 
2). It decreases from 200 Hz to 50 Hz if e = 10"'' and 
from 700-1000 Hz to 200 Hz if £ = 10"^. Consequently 
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Fig. 1. frequency distribution of detectable pulsars if £ = 
10-5 
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Fig. 2. frequency distribution of detectable pulsars if £ = 
10-6 

we could keep 90% of the maximmn detection probability 
while exploring less than 7.5% of the frequency space with 
an advanced detector versus 25% with interferometers of 
the first generation. Another important result concerns 
the lower bound of the detector frequency band. At the 
begining, this value doesn't play much part but the im- 
provement of the sensitivity would step by step allow to 
recover the all "potentially detectable" population con- 
centrated at the lowest frequencies. 

As already mentionned one of the major complication 
when searching for pulsars over long period of time arises 
from the intrinsic spindown which spreads power across 
many frequency bins. A model of the intrinsic frequency 
evolution is thus needed to remove these effects before per- 
forming the Fourier transform. Estimates of the spindown 
parameters can be derived from observations of radio pul- 
sars (Brady et al., 1998 and 2000 for instance) but strong 
selection effects make them non representative of the de- 
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Table 2. parameters of the detectable frequency distribu- 
tion for e = 10~^ (left value) and e = 10~^(right value) for 
interferometers of the first generation (LIGO and VIRGO) 
and for an advanced interferometer: the maximum, the 
width and the frequency band around the maximum where 
we find 90% of detection 



max (Hz) width (Hz) 90% band (Hz) 

VIRGO 180 - 480 200 - 700 [0-1200]- [0-2500] 

LIGO I 135 - 435 200 - 1000 [0-1000]- [0-2600] 

Advanced 45 - 85 50 - 200 [0-300]- [0-2500] 



Table 3. parameters of the period derivative distribution 
for £ = 10~^ (left value) and e = 10~^(in parenthesis) for 
interferometers of the first generation and for an advanced 
interferometer; the maximum, the 0.9 and 0.99 quantiles 
(lines 3 & 5) and the and the corresponding period shifts 
in nanosecond after one year of observation (lines 4 & 6) 



LIGO/VIRGO Advanced 



max -17.5 (-17.2) -16.8 (-16.4) 

logPg„% -16.7 (-16.4) -16.0 (-15.6) 

AP9o% (ns) 0.63 (1.26) 3.15 (7.93) 

logPgg% -16.1 (-15.8) -15.4 (-15) 

AP9g%(ns) 2.51 (5.00) 12.65 (31.56) 



T 




X = log Pdot 

Fig. 3. period derivative distribution of the detectable 
pulsar population if e = 10~^ for the sensitivity of in- 
terferometers of the first generation (plain line) as well as 
for an advanced interferometer (dashed line). The dot line 
represents the expected distribution of the 'potentially' 
detectable population 

tectable population. In the following, we investigate the 
statistical properties of the spindown rate of detectable 
pulsars. The period derivative distribution of detectable 
pulsars as a function of the antenna sensitivity and for an 
ellipticity e = 10~® is shown in Fig. 3 and compared with 
the expected distribution of the 'potentially detectable' 
population. As there is no significant differences between 
the distributions of the first generation interferometers, 
only one plot is reported for both LIGO I and VIRGO. 

The period derivative probability function can be rep- 
resented quite well by a gaussian distribution with stan- 
dard deviation a — 0.6. The maximum (or the mean 
value) is around log(P) = —17.5 for interferometers of 
the first generation and increases to log(P) = —16.8 for 
an advanced interferometer. As the sensitivity improves, 
the distribution shifts toward the expected distribution of 
the 'potentially detectable' population which maximum 
is log(P) = —15.6. Changing the mean ellipticity from 
£ = 10~^ to £ = 10~^ acts the same way (see Table 3) 
The general belief is that detectable pulsars will be the 



youngest at the beginning of the operations but it is shown 
here that it is not necessarily the case. Actually, detectable 
pulsars are in majority the 'most stable'. That is very un- 
derstandable considering that the less is its spindown, the 
longer a pulsar remain in the high frequency range, which 
increase its probability to be detected. 

It is relevant to define the critical spindown rate ^90% 
below which 90% of detectable pulsars can be found. This 
parameter may be used to determine the number of bins 
that should be investigated in the frequency space to cor- 
rect spindown effects. Table 3 gives ^90% as well as a 
more conservative ^99%, for interferometers of the first 
and second generations and for ellipticities e = 10^^ and 
10~^. Columns AP9o% and APgc,% give an estimation of 
the maximal shift one should expect (at 90% and 99% of 
confidence) after one year of observation (equation 7 with 
T = 1 yr) 

AP9o% = P90% X T (7) 



2.4. Conclusions 

In this paper, the distribution in the frequency space, of 
pulsars with a gravitational strain high enough to be de- 
tected by the future generations of laser beam interfer- 
ometers was investigated. Until detectors become able to 
recover the entire population, the 'detectable' frequency 
distribution will be very dependent on the detector sen- 
sitivity. If £ = 10~^, the optimal frequency is around 450 
Hz for interferometers of the first generation and shifts 
toward 85 Hz for an advanced detector. An encouraging 
result for future detection stategies is the significant nar- 
rowing of the optimal band of frequency which follows the 
improvement of the sensitivity: with an advanced detector 
we could keep 90% of the maximum detection probabil- 
ity while exploring less than 20% of the frequency space 
(7.5% £ = 10~^). In addition, we showed that in most 
cases 'detectable' pulsars have very stable rotations, mak- 
ing them easier to track in the frequency space. The crit- 
ical spindown rate below which 90% of detectable pul- 
sars can be found is of the order of a few nanosecond per 



Tania Regimbau and Jose Antonio de Freitas Pacheco: Searching GW from pulsars.. 



5 



year for interferometers of the first generation (and a tens 
of nanosecond per year for an advanced interferometer). 
These informations, integrated into detection algorithms 
may be usefuU to choose the optimal way to start blind 
searches, increasing our chances to detect these objects. 
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